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We study the Schrodinger functional coupling for lattice Yang-Mills theory coupled to an improved bosonic 
spinor field, which corresponds to QCD with minus two light flavors. This theory serves as a less costly testcase 
than QCD for the scaling of the coupling. 



1. Introduction 

The non-perturbative computation of the 
strong coupling constant as from hadronic scales 
to perturbative high energy scales is one of the 
main interests of the ALPHA collaboration. For 
this computation, a finite volume renormaliza- 
tion scheme is used in which the coupling runs 
with the space-time volume ||^. The scale evolu- 
tion of the coupling is computed recursively with 
the help of the step scaling function, which is an 
integrated variant of the beta function for finite 
changes of the scale. Furthermore, Schrodinger 
functional boundary conditions and 0(a) im- 
provement are employed. An important property 
of the step scaling function is that it can be calcu- 
lated by simulating pairs of lattices with different 
sizes L for the same a and taking the continuum 
limit. 

After this programme had been successfully 
used in the quenched approximation first re- 
sults for full QCD with two flavors are now avail- 
able Unfortunately, full QCD is notoriously 
costly to simulate, and thus the data for the step 
scaling function do not reach very close to the 
continuum limit. In order to learn more about 
the extrapolation to the continuum and about the 
effect of 0(a) improvement, we have studied the 
step scaling function in a Yang-Mills theory cou- 
pled to a bosonic spinor field, which has a local 
interaction and is therefore much cheaper to sim- 
ulate. It corresponds to setting the number of 
flavors iVf — —2 in the partition function, and 
is therefore closely related to full QCD in per- 
turbation theory. In the literature, this model is 



known as bermion model y . Our focus here how- 
ever is not on extrapolating results from negative 
to positive flavor number. 

2. Bermion model 

For a detailed discussion of the lattice setup, 
boundary conditions and the algorithm we have 
used, we refer to Q and references therein. 

The Schrodinger functional, as the partition 
function of the system, is an integral over all 
gauge and quark fields which fulfill the given 
boundary conditions. After integrating out the 
quark fields, it is 



Z = e 



-r 



D[U]e~'^<' det(MUf)^f/2 (^^^ 



with a gauge action Sg. For N{ = —2, the deter- 
minant can be written as 

Z= I D[U]D[(P+]D[4>]e~'^^~'^'', 



(2) 



with a now local bosonic action 
5b[[/,0]=a4V|(M</>)(x)r 



(3) 



*Talk given by B. Gehrmann at Lattice 2001, Berlin 



The fields 0(a;) carry color and Dirac indices. 

We have chosen the bulk improvement coeffi- 
cient Csw by extrapolating the non-perturbative 
data for A^f = § and A^f = 2 § in A^f. An 
explicit calculation of Csw along the lines of these 
references for the most critical parameters in our 
simulations has proven a good accuracy of this 
extrapolation. The boundary improvement coef- 
ficients Ct and Cf have been set to their perturba- 
tive values. 
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Figure 1. Cost for the measurement of (t(0.9793) 
for Wilson and improved bermions. 

As the improved bosonic action depends 
quadratically on each hnk Uxfi through the clover 
term, it is difficult to employ finite step size algo- 
rithms such as the combination of heatbath and 
overrelaxation steps used in [||. Because of the 
large additional cost of hybrid Monte Carlo algo- 
rithms, we have decided to use global heatbath 
steps for the bosonic fields and local overrelax- 
ation steps with respect to the unimproved action 
for the gauge fields. The clover term is taken into 
account by an acceptance step. The acceptance 
rate has turned out to be high enough (about 
70 %) for such an approach. Nevertheless, there 
is a significant overhead due to the calculation of 
the action difference necessary for the acceptance 
step, resulting in a cost factor of about 12 com- 
pared to unimproved Wilson bermions. On the 
other hand, comparing with data from im- 
proved bermions are about a factor of 10 cheaper 
than dynamical fermions, with a better scaling in 
a. In figure H we see a scaling of the cost with 

3. Results 



L/a = 16 at u = 1.5145 we have also used a 
crate of APEmillc. 

Figure § shows the results for the step scaling 
function plotted against (a/L)^. Within the er- 
ror bars, no linear dependence an a/L is visible, 
and therefore an extrapolation linear in (a/ L)'^ is 
justified. The extrapolated values are consistent 
with perturbation theory, and their error bars are 
of the same size as the perturbative 3-loop con- 
tribution. 
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We have computed the step scaling function 
T,{u,a/L) for the couplings u — = 0.9793 
and u = 1.5145 and lattice sizes L = 4,5,6,8. 
Most simulations were done on APEIOO parallel 
computers with up to 256 nodes with 50 MFlops 
each. For the simulation of the largest lattice 



Figure 2. Step scaling function for improved 
bermions for the couplings u — 0.9793 and u = 
1.5145 with fits linear in [a/L)"^. Also shown is 
the extrapolated continuum value and the 2- and 
3-loop results. 



3 



C\3 

I 

II 



Improved and Wilson bermions 




jA-loop PT . 1 








K 2-loop PT 









a/L 

Figure 3. Results for the step scaling function at 
u = 0.9793, together with a quadratic fit under 
the constraint of universahty. 

In the step scaling function was already 
computed in the unimproved bermion theory for 
u — 0.9793. Figure H shows these results together 
with the data after implementing improvement. 
The linear cutoff effects for this observable are of 
the order of a few percent. In the plot, both data 
sets are fitted under the constraint of universal- 
ity, i.e. that their continuum limit agrees. This 
fit is linear plus quadratic in a/L for the Wilson 
bermion data and quadratic in a/L in the im- 
proved data. Although the additional input from 
the Wilson data is included, the joint continuum 
limit crcombinod(0.9793) = 1.1059(43) is almost the 
same as the value cTimprovcd (0.9793) = 1.1063(46). 
A linear plus quadratic fit in a/L of the unim- 
proved data alone would have given the contin- 
uum result Crunimprovcd(0.9793) — 1.103(12). 

This indicates a success of the improvement 
programme. The main contribution to the cost 
for the calculation of a{u) comes from the largest 
lattice. When using an improved action, the lat- 
tice size needed for a reliable extrapolation to the 
continuum is smaller than without improvement. 
Even in the model used here, where the algorith- 
mic implementation implies a large overhead for 
the inclusion of the clover term, this leads to the 
improved case being more cost effective. In sim- 



ulations with dynamical fermions and algorithms 
like Hybrid Monte Carlo, the implementation of a 
clover term is possible with much lower overhead, 
and the advantage should therefore be higher. 
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